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Electron localization and polaron mobility in oxygen deficient as well as Li doped monoclinic 
tungsten trioxide have been studied. We show that small polarons formed in the presence of oxygen 
vacancy prefer the bipolaronic configuration whereas the configuration is 

found to be metastable. Our calculations suggest that bipolarons are tightly bound by the vacancy 
and therefore largely immobile. On the contrary, polarons formed as a result of Li intercalation can 
be mobile, the activation energy for polaron jumping in this case varies between 98 and 124 meV 
depending on the crystallographic direction. The formation of bipolarons in Li — WOs 

is possible. When situated along [001] the bipolaronic configuration is 8 meV lower in energy than 
two separate polarons. 


I. I. INTRODUCTION 

Phenomenon of electrochromism is found in many ma¬ 
terials such as: organic compounds [T], polymer compos¬ 
ites ElEl and semiconductors [4] . One of the most stud¬ 
ied materials of this class is tungsten trioxide due to a for¬ 
tunate combination of its corrosion resistivity, sufficient 
light absorption and satisfactory mechanical properties 
that make it an attractive material for advanced applica¬ 
tions [ME] such as ’’smart windows” technology mm- 

The first mentioning of WOs coloration upon ion inser¬ 
tion dates back to the 19th century m- In the 20th cen¬ 
tury the interest to this material and its electrochromic 
properties revived and resulted in a large number of stud¬ 
ies inHisi- It is well known by now that ’’Prussian blue” 
coloration in the crystalline and amorphous phases of 
tungsten trioxide appears upon hydrogen or alkali metal 
(Li, A^a) intercalation l9l ITOi ITBllTS] . 

Tungsten trioxide is a wide band-gap semiconductor 
(2.6-3.0 eV [171 [191 HO]) with perovskite like structure, see 
Eig.l. The oxide shows rich structural polymorphism 
mHH], whereby oxygen octahedrons surrounding tung¬ 
sten atoms vary their orientations and shape depending 
both on temperature and pressure. Usually, at room tem¬ 
perature both 5-WO^ (triclinic PI 233-290 K) and 7- 
IUO 3 (monoclinic P 2 \jn 290-350 K) phases coexist [T7|. 
Depending on the method of preparation samples can be 
characterized by various 7/J ratios, micro-structure and 
defectiveness. Under certain conditions a more exotic 
hexagonal h-WOs phase j^S] HB] or amorphous a-LEOs 
m phase can form. The structural variety of WO 3 of¬ 
fers a range of optical characteristics that is important 
for applications relying on light absorption [27H29]. 

Alkali metal intercalation can also cause structural 
transformation of IUO 3 . Eor example, in the case of 
LixWO^ several phase transitions from monoclinic to cu¬ 
bic symmetry with increasing Li content have been re¬ 
ported [30H32]. It has also been shown that for Li+ 
intercalation the degree of coloration increases propor¬ 
tionally to the lithium concentration whereas in the case 


of hydrogen insertion a saturation level of coloration is 
reached for a certain concentration [33]. Based on 
neutron diffraction data [HU (35) this difference between 
H and Li has been explained by the formation of O — H 
bonds upon hydrogen insertion. At the same time alkali 
metal ions show no tendency of binding to the oxygen 
atoms of the oxide matrix. 

Substoichiometric WO 3 -X also shows multifarious 
electrochromic behavior, demonstrating different degrees 
of coloration depending on the concentration of oxygen 
vacancies IMlISZj- Thus, experiments clearly indicate 
that both doping (P, Li) and oxygen vacancies can influ¬ 
ence WO 3 optical properties [33l [37H39] . 

Two general concepts of the colouration mechanism in 
tungsten oxide have been suggested: inter-valence charge 
transfer [4] and polaron models gniiii], the latter be¬ 
ing more commonly accepted. According to the polaron 
model the coloration is the result of photon absorption 
due to photon activated polaron hopping. A polaron 
can form in WO 3 due to electron localization at the d- 
states of tungsten, thus, turning into or . 
Light absorption can lead to the activation of the fol- 

lowing transitions; W®+i + W^^le 2 ^ 

+ w!ite 2 ^ + W,te 2 - We notice that more 

complex configurations of localized charge can form, for 
example, two bound polarons at neighboring sites, 1 U^+- 
1U^+, also called bipolarons (see Eig.lb,c for illustration). 

The presence of different localized charge states in 
WO 3 has been confirmed by numerous experiments 
ESj. The formation of poiarons in this oxide has aiso been 
studied theoreticaiiy by semi-empiricai modeis [501 EH 
and more recentiy by ab-initio methods empioying hy¬ 
brid functionais [53] . More studies of the eiectronic and 
structurai properties of perfect and defective WO 3 can 
be found in Refs. jsaisiHsii- Despite intensive research 
many questions regarding the properties of poiarons in 
WO 3 remain unciear. In this paper we address some of 
them, in particuiar, we study the formation and mobiiity 
of poiarons in two cases: 1 ) in the presence of an oxygen 
vacancy and 2) in the presence of Li. We investigate the 
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Figure 1. The structure of monoclinic 7 -VFO 3 phase without defect (a), with an oxygen vacancy (b), with two Li atoms c). 
Charge distributions are shown for the polaronic states. 


relative stability of polarons and bipolarons and calculate 
the energy barriers for polaron propagation in WO 3 . 


II. II. DETAILS OF CALCULATIONS 

The presented calculations were performed in the 
framework of DFT using the projector augmented-wave 
method, as implemented in VASP [58l |59]. The 
DFT + U [60] approach using generalized gradient ap¬ 
proximation (GGA) in the Perdew, Burke and Ernzer- 
hof {PBE) [H] parameterization was used to account 
for the exchange-correlation interaction. The U values 
{Uop = and Uwd = 6 eU) were determined using 
the linear response procedure implemented in Quantum 
ESPRESSO package (62] [63]. Hereafter we refer to this 
set of calculations as to DTF + Uw^Uq- We also per¬ 
formed the DFT + U calculations with Hubbard — U 
applied only to the d-states of W{DFT F Uw)- Ad¬ 
ditionally, we used the screened hybrid functional of 


Heyd, Scuseria, and Ernzerhof {HSE06) [^11115] with the 
Coulomb potential separation parameter, cj, of 0.2 and 
the HF mixing constant of 0.25. 

Our DFT -f U calculations were performed for the 
2 x 2 x 2 unit cell (192 O and 64 W) and 2 x 1 x 1 unit cell 
(48 O and 16 LF), both using a P-point centered 2 x 2 x 2 
k-point grid for the Brillouin zone integration. The 
2 x 1 x 1 unit cell was also used in the hybrid calcula¬ 
tions. The following states were treated as valence states: 
5p^, 5d^, 65 ^(LF) and 2s^,3p^(0) and ls^, 2 p^{Li). All 
the calculations were spin-polarized. The cut-off energy 
was 650 eV. 

The structure relaxation was done using the following 
scheme. Eirst, we setup an initial supercell using the ex¬ 
perimental parameters [23] then we relaxed the structure 
in two steps: 1 ) first we fixed the atomic positions and 
relaxed only the cell parameters; 2 ) next we fixed the 
cell volume and relaxed the atomic positions. The latter 
were relaxed until the Hellmann-Eeynman forces acting 
on atoms became smaller than 10“^ eV/A. We repeated 
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Table L The lattice parameters of the monoclinic cell obtained in the HSE06 and DFT + Uw^ Uo calculations together with 
experimental values |23) . In calculations angle /3 was set to the experimental value 90.89° |23| . 


Lattice parameters, A 

Experiment |23| 

HSE06 

{/(Op) = 9eV, U{Wd) = 6eV 

a 

7.30 

7.34 

7.41 

b 

7.54 

7.58 

7.65 

c 

7.69 

7.73 

7.81 


this two step relaxation routine a few times to ensure a 
stable result. We notice that during the relaxation we 
kept the symmetry of the lattice monoclinic. 

Here we study only uncharged oxygen vacancies, for 
which the vacancy formation energy, Ey^ can be calcu¬ 
lated as 


Ev = Ew^Osa, - Ew^Osx-i ~ 1/2-E’o2 5 (1) 

The barriers for polaron transitions were calculated us¬ 
ing both the linear interpolation scheme {LIS) as in our 
previous work [66] and the nudged elastic band method 
{NEB) [67|. NEB is an efficient method for studying 
atomic and molecular diffusion, phase transitions, chem¬ 
ical reactions as well as polaronic transitions in oxides 
[67HZI]. We used 10 intermediate configurations in LIS 
and 5 images for NEB. 


10 

5 

0 


-- -5 



§3 10 
B 5 

on 

‘S 0 
^ -5 

1 -10 

(D 

^ 10 
5 
0 
-5 
-10 


1 1 

1 1 1 1 

U ^=8 eV 

1 _i 

i_1_1_1 


1 1 1 

▲ 

1 1 1 1 

Uw =6 eV 

Uo =9 eV 

1 _ 1 _i 

i_1_1_1 


-1-1-1-1-1-1-j 

!-'-'-'- 

HSE 06 

_1_1_1_1_1_1_ 1 

i_1_1_1_ 


-3.5-3.0-2.5-2.0-1.5-I.O-0.5 0 0.5 I.O 1.5 2.0 


Energy relative to the E^, (eV) 


III. III. RESULTS AND DISCUSSION 

A. A. Electron localization in oxygen deficient 

'j — WO 3 

To properly model small polarons one needs a method 
adequately describing the on-site charge localization and 
lattice distortions around this site [saiziHzs]- Electron 
localization in correlated d— and /—metal oxides usually 
cannot be correctly described by standard DFT approxi¬ 
mations, which due to the presence of the self-interaction 
term tend to overestimate the stability of the delocalized 
solution. The Coulomb interactions can be corrected, for 
example, by means of the Hubbard—U parameter (75]. It 
can be treated as an adjustable parameter to better de¬ 
scribe some experimental data m or it can be obtained 
from restricted ab initio calculations m- Traditionally 
in oxides Hubbard — U is applied to the d— or /—states 
of cations, in some cases, however, the description of 
oxide property can be further improved by applying U 
also to the p—states of oxygen IZMll- Another way to 
cure the self-interaction problem of standard DFT is to 
employ the so-called hybrid functionals, where a certain 
amount of nonlocal Eock exchange is added [80|. Both 
approaches have successfully been used to describe elec¬ 
tron localization or small polaron formation in different 
oxides 123 is [71HZS1 mi- Hybrid functionals are more 
computationally demanding than DFT-\-U and therefore 
their usage is limited to rather small unit cells. 


Figure 2. DOSs calculated for the 2x1x1 supercell with 
an oxygen vacancy and associated bipolaron(see 

Fig.lb). The top, middle and bottom panels present DOSs 
obtained with DFT + Uw^ DFT -h Uw, Uo^ HSE06, respec¬ 
tively. The total DOSs are shown in grey, the polaronic Wd 
states are in blue. 


To find the optimal functional for the description of 
polarons in WOs we have performed three sets of cal¬ 
culations using the 2x1x1 unit cell (64 atoms), namely, 
DFT + Uvp, DFT + Up/, Uq and HSE06 (see Sec. H for 
details). As a test case we considered an oxygen vacancy 
with two unpaired electrons localized on both sides of the 
vacancy in the configuration (Fig.lb). 

First, within the DFT Uw approach we varied Uw 
from 2 to 12 eV in the manner described in Ref. |66| . We 
analyzed the localization patterns for different values and 
found that for Ueff = SeV and higher almost precisely 
one electron was localized at the d-orbitals of each of the 
two tungsten atoms. The local magnetic moment 

{fiw) for each W in this case was 0.96 iiB. Next, we per¬ 
formed calculations using DFT T Uw^Uq- In this case 
the local magnetic moment due to the electron localiza¬ 
tion at the H/^-states was 0.86 iiB. Finally, our HSE06 
calculation resulted in fiw = 0.64 laB. 

The density of states {DOSs) obtained in the three 
types of calculations are compared in Fig.2. The band 
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Table II. Bond distances and angles obtained for the monoclinic WO3 structure using DFT + Uw,Uo- The corresponding 
experimental values are shown in parentheses ( |23|b The oxygen position notations are as shown in Fig. la. The subscripts 
x,y,z are added for convenience to indicate the orientation of the corresponding pair in the crystal. 


Oxygen position 


WO-bondl, 0 

WO-bond2, () 

ou, [100] 

160.1, (154.8) 

1.90, (1.89) 

1.89, (1.81) 

02 ,, [100] 

160.7, (155.7) 

1.90, (1.95) 

1.89, (1.91) 

03y, [010] 

163.6, (161.5) 

1.84, (1.82) 

2.03, (2.02) 

04y, [010] 

163.4, (160.1) 

2.04, (2.06) 

1.84, (1.78) 

05,, [001] 

170.5, (164.2) 

2 .10, (2.16) 

1.80, (1.71) 

06,, [001] 

169.7, (157.3) 

1.81, (1.81) 

2.09, (2.11) 


gaps are expectedly underestimated by both versions of 
the DFT + U calculations: Eg {DFT + [/vp)=l-43 eV, 
Eg{DFT + Uw^ Uo)= 1.64 eV that should be compared 
to the experimental value of 2.8 eV mi. At the same 
time, HSE06 overestimates the gap: Eg{HSE0Q) = 3.19 
eV. In the case of DET Uw the localized polaronic d- 
states are situated at the upper edge of the valence band 
whereas for HSE06 and DET FUw^Uo these states are, 
respectively, 0.8 eV and 2.5 eV above the valence band 
top. The position of the localized states with respect to 
the conduction band edge is similar for both approxima¬ 
tions. We notice that the position of the polaronic peak 
in the gap obtained by HSE06 is in good agreement with 
previously reported B3LYP results m- 

The magnetic state of the obtained bipolaron is de¬ 
generate according to HSE06 whereas in the case of 
DET + Uw^Uo the antiferromagnetic coupling of local 
spins is favored by about 20 meV over the ferromagnetic 
one. Using DET -b Uw^Uo we also examined DOSs 
calculated for Li — WO^ using the 2 x 2 x 2 unit cell and 
found them to be very similar to those calculated for the 
vacancy case. The only noticeable difference was a 40 
meV shift of the polaronic peak down in energy. 

The lattice parameters obtained using DET + Uw^ Uq 
and HSE06 together with experimental data are shown 
in Table I. They demonstrate reasonable agreement with 
the HSE06 values being a bit closer to the experimental 
parameters. Based on the analysis of DOSs^ localization 
behavior and structures obtained in the three sets of cal¬ 
culations we conclude that the DET FUw^Uo approach 
is the most appropriate for studying polarons in WOs 
as it shows good agreement with the results of HSE06 
and, at the same time, allows us to treat large supercells. 
Therefore, further on we study polaron mobility in WOs 
using the DET FUw^Uo approach. 


B. B. Analysis of oxygen vacancy positions and 
vacancy trapped polarons. 

The monoclinic phase of WOs (2ci/n) has two tung¬ 
sten and six oxygen non-equivalent atoms [23], as il¬ 
lustrated in Fig.la. The structure is rather complex, 
each O atom has two neighboring W atoms but the 
W — O — W chains differ in all the directions. Along 


[100] we find ... -01-W1-02-W2-01- ... and 
... — 02 — lUl — 01 — 1U2 — 02 — ... chains, along [001] 
... - 06 - VFl - 05 - VF2 - 06 - ... and ... -05-VFl- 
06 — 1U2 — 05 — ... chains interconnected in the [010] 
direction by alternating ... VFl — 03 — Wl — 03 — VFl — ... 
and ... — kF2 — 04 — kF2 — 04 — kF2 — ... chains (for 
notations see Fig.la). The W — O bonds in all the direc¬ 
tions have alternating shorter and longer distances and 
slightly different O — VF — O angles, which are shown in 
Table II, both calculated with DET FUw^Uo and from 
the experimental data [23]. 



Figure 3. Vacancy formation energies for different vacancy 
positions and two electronic configurations: and 

pp5+ _ ^ 5 + 

energies are shown with respect to the 
lowest formation energy, the one for the 05^ position, with 
the energy of 2.36 eV. See also Table III. 


Using the structure obtained in DET FUw^Uo we 
calculated the vacancy formation energetics for the six 
oxygen positions (see Fig.la) and two bipolaronic config¬ 
urations: kF^+ — kF^+ and kF^+ — kF^+ bipolaron (see 
Fig.lb). Both electronic configurations could be stabi¬ 
lized for all the vacancies. The results are presented in 
Fig.3 and Table III. The lowest vacancy formation energy, 
2.36 eV, is found for the oxygen vacancy at site 05z with 
bipolaronic configuration — and the highest en¬ 
ergy, 2.90 eV, for the 01 position in the con- 
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Table III. Vacancy formation energies, Ey, and relevant W-W distances, R, for non-equivalent oxygen sites and — (see 
Fig.lc) and configurations. The numbers in parentheses show the corresponding W — W distances in the perfect 

WO 3 cell. 


Oxygen crystallo¬ 

graphic positions 

E4eV),W^+ - W^+ 

E4eV),W^+ - VC^+ 

Ri),W^+ -W^+ 

R{),W^+ - W^+ 

ou, 

2.79 

2.90 

3.94, (3.72) 

3.87 

02a:, 

2.77 

2.86 

3.92, (3.74) 

3.86 

03y, 

2.66 

2.78 

4.01, (3.83) 

3.90 

O^yi 

2.62 

2.73 

4.05, (3.84) 

3.92 

05., 

2.36 

2.51 

4 . 35 , ( 3 . 91 ) 

4.19 

06., 

2.40 

2.54 

4.31, (3.89) 

4.14 



Figure 4. Energies along the transition path between the 
and configurations for three vacan¬ 

cies {02x, 04y and 05z) as estimated by NEB. The energies 
are given with respect to that of in the cor¬ 

responding direction. For comparison for 05z we show the 
results obtained by the linear interpolation method {LIS) 
|66| . Charge distributions corresponding to the initial x=0 
(pp5+ _ ^nal x=l (VF^^) transition points are 

shown for illustration. 


figuration. The formation energies for the 
configurations are systematically higher than those for 
the ones (Table III). The vacancy formation 

energy, Ey^ depends on the orientation of the bipolaron in 
the lattice, Ey averaged for each direction can be sorted 
in the following order: F^^[100] > F^^[010] > F^^[001]. 
These results are in agreement with previous studies car¬ 
ried out by standard DFT EZl, where the oxygen atoms 
with the W — 0 bonds along the [001] direction were indi¬ 
cated as favorable sites for vacancy formation. Using the 
B3LYP hybrid functional Wang et al. [53] also reported 
the formation of the configuration aligned 

with [001]. 

The data in Table III show that the vacancy forma¬ 
tion energy decreases as the W — W distances increase 
and lowest Ey corresponds to the largest W — W separa¬ 
tion. The formation of bipolarons {W^~^ — W^~^) leads to 


a noticeable increase in the W — W bond lengths due 
to Coulomb repulsion. These local lattice distortions 
are, however, different for different crystallographic di¬ 
rections. In the [100] and [010] directions the W — W dis¬ 
tance increase is of about 0.20 A whereas for the 06z and 
05z vacancy positions ([ 001 ]) it is 0.42 A and 0.44 A, re¬ 
spectively. We see a similar trend for the po- 

laronic configurations (Table III), however, the changes 
of the W — W distance are smaller in this case. These 
results indicate certain softening of phonon modes along 
[001]. They also indicate that in the case of monoclinic 
W O 3 the electronic contribution to the energy of vacancy 
formation dominates over the elastic contribution. 

Our calculations suggest that the bipolaronic states 
associated with an oxygen vacancy are tightly trapped 
as all our attempts to localize charge away from the va¬ 
cancy were unsuccessful. Therefore, we can expect that 
in the presence of oxygen vacancies bipolarons are im¬ 
mobilized in the lattice, however, the transitions between 
vacancy bound 1 U^+ - 1 U^+ and 1 U^+ - 1 U^+ could still 
contribute to the coloration effect in IUO 3 . To study 
the transition between these two configurations we chose 
oxygen vacancies with lowest Ey in each crystallographic 
direction, namely, 02^,, OAy and 05z> Our results, pre¬ 
sented in Fig.4, show that the transition barriers vary 
from 150 to 230 meV, being the smallest for the [100] di¬ 
rection and largest for [001]. This result is in agreement 
with similar trends shown by the energy differences be¬ 
tween the and configurations: 

02x, (90 meV) < OAy, (110 meV) < 05z, (150 meV) 
and the W — W bond length differences: 02^, (0.06 ) < 
04^, (0.13 ) < 05;^, (0.16 ). We also notice that the 
NEB relaxation results in lower barriers compared to 
those obtained with LIS (Fig.4). This is understand¬ 
able as NEB allows one to optimize structures along the 
transition path whereas LIS uses static configurations. 

Lattice distortions around polarons and especially 
bipolarons can be quite substantial reaching far beyond 
first coordination shells. In the case studied here we find 
that the local distortions around the bipolaron —vacancy 
complex spreads as far as 10-11 A in the direction of the 
alignment. The 2x2x2 cell used in our 
calculations is big enough to accommodate such a defect. 
Our results, however, highlight the importance of using 
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large unit cells in modeling polarons and bipolarons at 
least if one wishes to study isolated quasiparticles of the 
sort. 


C. C. Polaron formation and transitions in Li 
doped monoclinic 7 phase. 

Polarons can also form in the presence of impurities, 
dopants or due to electron injection. Here we report our 
results of the polaron formation and mobility in Li doped 
monoclinic H^Os. First we consider one lithium atom 
situated in an octapore in the supercell of 256 atoms, 
i.e. VFOsI/io.oie- This is a realistic concentration often 
studied in experiment [82] . After the relaxation Li stayed 
close to the octapore center (Fig Ic), no bond formation 
with the surrounding oxygen atoms was detected. It was 
possible to localize polaron at different tungsten 

sites. The association energy of Li~^ and is of the 

order of 30 meV, which was estimated by calculating the 
energies of separated by 3.92 A and 7.71 A. 

Further, to study polaron mobility we used tungsten 
sites away from Li 6-8 A) simulating a free polaron 
propagation in the WO^ matrix. Fig.5 shows the energy 
barriers for polaron jumping between two neighboring W 
sites in the three crystallographic directions. As expected 
the LIS barriers are again higher than those obtained 
using NEB. The NEB barriers vary from 98 for [001] 
to 124 meV for [100]. These values are in surprisingly 
close agreement with the activation energy for polaron 
hopping estimated from basic models using experimental 
information, 180 meV [82] . 

If the concentration of lithium ions and polarons in¬ 
creases the number of possible configurations of local¬ 
ized charges becomes larger. We study the possibility 
of a bipolaron formation by considering two Li atoms 
in the 256 atom IFO 3 supercell (lF 03 l/io.o 32 )- Next to 
the bipolaronic configurations we also con¬ 
sider the formation of two separate polarons and 

one . In Fig .6 we show the energies for different 
electronic configurations and transitions between them 
obtained using NEB. We start from the situation when 
both electrons donated by the two Li atoms are delo¬ 
calized in the oxide matrix (x=0. Fig. 6 ), then we local¬ 
ize one polaron that happens without any barrier 

(x= 2 . Fig. 6 ). After that we localize the second electron 
at the next nearest (NN) W site from first which 

neither involves any barrier(x=3. Fig. 6 ). The energies 
of such NN configurations are very close to each other 
for all the orientations (Fig. 6 ). Next we model a tran¬ 
sition of one from the NN site to the nearest (N) 

W site, thus, forming a bipolaron along the [100], [010] 
and [001] directions (x=3 ^ x=5. Fig. 6 ). This transition 
involves a barrier of about 90 meV. The resulting bipo¬ 
laron situated along [ 001 ] has the lowest energy not only 
among the bipolarons but among all the considered con¬ 
figurations (Fig. 6 ). The [001] bipolaron is 8 meV lower 
in energy than the [ 001 ] NN configuration. At the same 


time, the energies of bipolarons oriented along [ 100 ] and 
[010] are about 20 meV higher than those of the NN 
configurations. We notice that the localization of two 
electrons at one tungsten site is energetically un¬ 

favorable as this state lies more than 300 meV above any 
pp 5 +_pj 75 + configuration. Therefore, our results suggest 
that in Li doped WOs one can expect the formation of 
polarons and bipolarons whereas the 

presence of the state, often detected in experiments 
[29l |33l |36l Sni SHI S9], most likely can be attributed to 
oxygen vacancies, around which the state is found 

to be metastable. 



Figure 5. Energies along the transition paths of from 
one tungsten site to the other in the case of Li doped 7 phase. 
The energies are given with respect to the energy of . The 
solid lines show the NEB energies for all the three crystallo¬ 
graphic directions. For comparison, for [100] we also show the 
energies obtained in LIS calculations (dashed line). Charge 
distributions for at the initial x = 0 and final x — 1 

transition points are shown for illustration. 


IV. IV.CONCLUSION 

We have studied the applicability of iL5'F^06 and two 
DET + U approaches to modelling of polarons in IFO 3 . 
We have shown that the DET + U approach with two 
Hubbard — U parameters determined and applied simul¬ 
taneously to the d-W {Uw=^ dV) and p-0 {Uo=9 eV) 
states provides a proper description of polaron forma¬ 
tion and transitions in IFO 3 . At the same time, this 
method allows us to calculate large unit cells necessary 
for adequate modelling of polarons and bipolarons. Us¬ 
ing this approach we have studied the vacancy energet¬ 
ics for the six non-equivalent oxygen positions and two 
electronic configurations: and . 

The bipolarons situated along [001] are the 

most favorable electronic configurations around vacan¬ 
cies. We have also studied polaron and bipolaron for¬ 
mation and energetics for Li doped IUO 3 . In this case 
the bipolaronic configuration aligned with 
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Figure 6. Energies along the transition paths calculated using NEB for the case of two intercalated Li atoms (Fig.lc). The 
energies are shown with respect to the lowest energy configuration situated along [001]. The following polaronic 

transitions are shown: from x=0 (two delocalized electrons) ^ x=2 (one electron localized as the other is delocalized) 

then ^ x=3 (two as next nearest neighbors) and further ^ x=5 — bipolaron configurations). These transitions 

are calculated along [100], [010] and [001]. The energy of the configuration is marked by the dashed line. Two dotted 

lines indicate the energy difference between the IF^^ — IF^^ bipolaron and NN configuration situated along [001] direction. 
Charge distributions for the IF^^ polarons at x=2, 3 and 5 are shown for illustration. 


[001] has again the lowest energy but winning only 8 meV 
over two separated 1F^+. Our results suggest that the 
polarons formed due to oxygen vacancies are immobile, 
at the same time, the 1F^+ state is metastable and the 
1^5+ _ ^5+ ]F^+ — 1F^+ transition is possible with a 

barrier of 150 meV. On the contrary, polarons formed in 
Li doped tungsten oxide are mobile with the minimum 
activation energy (98 meV) in the [001] direction. The 
1F^+ state is 300 meV higher in energy than any stud¬ 
ied 1F^+ configuration and, therefore, 1F^+ is unlikely to 
form in perfect Li — WOs without presence vacancies or 


similar structural defects. 
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